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ABSTRACT 

We present high-resolution continuum images of the W51e2 complex processed from archival data of the Sub- 
millimeter Array (SMA) at 0.85 and 1.3 mm and the Very Large Array at 7 and 13 mm. We also made line images 
and profiles of W51e2 for three hydrogen radio recombination lines (RRLs; H26a, H53a, and H66a) and ab- 
sorption of two molecular lines of HCN(4-3) and CO(2-l). At least four distinct continuum components have 
been detected in the 3" region of W5 le2 from the SMA continuum images at 0.85 and 1 .3 mm with resolutions of 
0.3" x 0.2" and 1 .4" x 0.7", respectively. The west component, W5 Ie2-W, coincides with the ultracompact HII re- 
gion reported from previous radio observations. The H26a line observation reveals an unresolved hyper-compact 
ionized core (< 0.06" or < 310 AU) with a high electron temperature of 1.2 x 10 4 K, with the corresponding 
emission measure EM> 7 x 10 10 pc cm" 6 and the electron density N e > 7 X 10 6 cm" 3 . The inferred Lyman con- 
tinuum flux implies that the HII region W51e2-W requires a newly formed massive star, an 08 star or a cluster of 
B-type stars, to maintain the ionization. W51e2-E, the brightest component at 0.85 mm, is located 0.9" east from 
the hyper-compact ionized core. It has a total mass of ^140M Q according to our spectral energy distribution 
analysis and a large infall rate of > 1.3 x 10~ 3 M Q yr" 1 inferred from the absorption of HCN. W51e2-E appears 
to be the accretion center in W51e2. Given the fact that no free-free emission and no RRLs have been detected, 
the massive core of W51e2-E appears to host one or more growing massive proto-stars. Located 2" northwest 
from W51e2-E, W51e2-NW is detected in the continuum emission at 0.85 and 1.3 mm. No continuum emission 
has been detected at A > 7 mm. Along with the maser activities previously observed, our analysis suggests that 
W5 le2-NW is at an earlier phase of star formation. W5 Ie2-N is located 2" north of W5 Ie2-E and has only been 
detected at 1.3 mm with a lower angular resolution (~ 1"), suggesting that it is a primordial, massive gas clump 
in the W51e2 complex. 

Subject headings: HII regions - ISM: individual objects (W51e2) - stars: formation 



1. INTRODUCTION 

Massive stars are formed in dense, massive molecular cores. 
Detailed physical processes in star-forming regions have not 
been well studied until recent high-resolution observations were 
available at submillimeter wavelengths. High angular resolu- 
tion observations are necessary to unveil the physical environs 
and activities of the individual sub-cores as well as their impact 
on the overall process of massive star formation. 

W51 is a well-known com plex of HII regio ns with about 
1° area in the Galactic plane dWesterhouli ri958). W 51A, also 



calle d G49.5-0.4, is the most luminous region in W5 1 (IKundu & Velu 
Il967l) . Several discrete components have been found from W5 1 A 
and re ferred alphabe t ically as a-i in th e right ascension (RA.) 
order (lMartinlfl97l lMehringerf[l99l . W51e is one of the 
brightest regions in W51A. High-resolution centimeter obser- 
vations have revealed t hat W51e consists of several ultracom- 
pact (UC) HII regions dScottHl978tlGaume et allll993l) . among 
which the UC HII region of W51e2 appears to be the bright- 
est. This region is severely obscured by the dust, and no IR 
detection has been made at wavelengths shorter than 20 



2004: lKetoet al.ll2008l) . 

Observations of hydrogen radio recombination lines (RRLs) 
were carried out toward W5 1e2 ( Mehring erl 1 1 994t iKeto et all 
120081: iKeto & Klaassenll2008l) . iMehringerl d!994h failed to de- 
tect the H92a line in this region, which might be attributed 
to the large optical depth and/or the large pressure broaden- 
ing of the UC HII region at 3.6 cm. Based on the observations 
of H66a, H53a, and H30a, IKeto et alj d2008l) argued that the 
pressure broadening is responsible for the broad-line widths ob- 
served in the RRLs. 

Ho & YoungH995lZhang & Hoi 



.Obser vations of molecular lines (e.L ^ - 

IZhang et alll99llYoung et alll998hlSollins et al.l l2004) 
showed evidence for infall (or accretion) gas within 5 " (< 0.2 pc) 



Genzel et al.ll982l) . The inferred distance to W5 le2 is ~5. 1 kpc 



dXu et al.l l2009): thus, 1" corresponds to a linear scale of about 



5100AU. 

Based on the radio continuum observations, spectral energy 
distribution (SED) analyses of W5 1 e2 suggest that the emission 
source from this region consists of two components, namely, (1) 
thermal emission from a cold d ust component and (2) free-free 
emission from a hot HII region (Rudolp h et al.ll99 0: SoHin s~et al.1 
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around W51e2. A possible rotation was suggested by Zhang & Hoi 
(1997) on the basis of fitting the position-velocity (PV) dia- 
gram of the NH3(3,3) absorption line at 13 mm. A spin-up ro- 
tation with an axis of position angle (PA.) «20° was further 
sugge sted based on the v elocity gradient of CH?CN(83-73) at 
2 mm dZhang et al.lll998l) . However, based on the velocity gra- 
dient observed from the H53a line. lKeto & Klaassenl d2008l) in- 
terpreted their results as evidence for rotational ionized accre- 
tion flow around the UC HII region and derived a rotational axis 
of P.A.Ri -30°. In addition, a possible bipolar outflow along 

the northwest and southeast (NW-SE) directions was s uggested 

based on the observations of the CO(2- 1 ) line emission dKeto & Klaa ssen 
2008), which appears to be perpendicular to the axis of the ro- 
tational ionized disk. However, the southwest (SW) elongation 
of the UC HII region in Gaume et al.'s (1993) observations was 
explained as a collimated ionized outflow. On the other hand, 
proper-motion measurements of H2O/OH masers appeared to 
favor a hypothetica l model with an outflow or expan ding gas 
bubble in W5 le2 dlmai et alj|2002l: iFish & Reidll2007h . 

Thus, improved images and comprehensive analysis of high- 
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Table 1 

Observation parameters for W5 1e2 



Parameters 



SMA: 0.85 mm 



SMA: 1.3 mm 



VLA:7mm VLA:13mm 



Observation date 
Array configuration 
Pointing center R.A.(J2000) 
Pointing center decl.(J2000) 
Frequency(GHz) 
Bandwidth 
On-source time (hr) 
System temperature(K) 
Bandpass calibrators 
Phase calibrators 

Flux calibrators 



2007 Jun. 18 

Very extended (8 ants) 

19:23:43.888 

+14:30:34.798 

343, 353 

2.0+2.0 (GHz) 

2.10 

150-500 

J1229+020, J1751+096 
J1733-130, J1743-038 
J2015+371 
Callisto 



2005 Sep. 1 
Extended (6 ants) 
19:23:43.895 
+14:30:34.798 
221, 231 
2.0+2.0 (GHz) 
4.75 

100~250 
3C 454.4 

J175 1+096, J2025+337 
Ceres 



2004 Feb. 14 2003 Sep. 29 

AB BC 

19:23:43.918 19:23:43.918 

+14:30:28.164 +14:30:28.164 

43 22 

12.5 (MHz) 12.5 (MHz) 

5.90 3.93 



3C 84 
J1923+210 

3C 286 



3C 84 

J 1923+2 10 

3C 286 




resolution observations at the wavelengths from radio to sub- 
millimeter are needed to reconcile the differences in the in- 
terpretation of the results from the previous observations of 
W51e2. In this paper, we show the high-resolution contin- 
uum images with the Submillimeter Array (SMA) at 0.85 and 
1.3 mm and with the Very Large Array (VLA) at 7 and 13 mm, 
as well as line images and profiles of the hydrogen recombina- 
tion lines (H26a, H53a, and H66a) and the molecular lines of 
HCN(4-3) and CO(2-l). In Section 2, we describe the details 
of the data processing and results. In Section 3, we model the 
individual components on the basis of the high-resolution ob- 
servations of W51e2 with the SMA and VLA. The nature and 
astrophysical processes of star formation cores in W51e2 are 
discussed. A summary and conclusions are given in Section 4. 

2. DATA REDUCTION AND RESULTS 

We processed the archival data from the SMA and VLA ob- 
servations of W51e2 and constructed both the continuum and 
spectral line images. Details in calibration and imaging are dis- 
cussed below along with the relevant parameters that are sum- 
marized in Table [TJ 

2. 1 . The SMA continuum images at 0.85 and 1.3 mm 

The observations of W5 le2 at 0.85 mm were carried out with 
the SMA in the very extended array (eight antennas) and at 
1.3 mm in the exten ded array (six a ntennas). The data were cal- 
ibrated in Miriad dSault et alJI 1995b following the reduction in- 
structions for SMA data0. The system temperature corrections 
were applied to the visibility data. The bandpass calibration 
was made with strong calibrators for each of the two data sets 
(TableQJi. The flux density scale was determined using Callisto 
at 0.85 mm and Ceres at 1.3 mm. 

Corrections to the complex gains at 1.3 mm were made by 
applying the solutions interpolated from the two nearby QSOs 
J1751+096 and J2025+337. For the high-resolution data at 
0.85 mm, the visibility data are initially calibrated with a model 
of W5 le2 derived from observations in a low angular resolution 
at 0.87 mm. The residual phase errors were further corrected 
using the self-calibration technique. The final continuum im- 
ages were constructed by combining all the line-free channels 
in both upper sideband and lower sideband data. 

Figure QJ shows the high-resolution (0.3" x 0.2") continuum 
image of W51e2 at 0.85 mm. The complex of W51e2 has been 
resolved into at least three bright, compact components (see 
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FIG. 1. — High-resolution 0.85 mm continuum image of W51e2 ob- 
served with the SMA. Contours are ±5a X 2" (n = 0, 1, 2, 3, ... and 
cr = 7.11 mJy beam" 1 )- The FWHM beam (0.31" X 0.22", P.A.=60.6°) is 
shown at bottom right. Three continuum emission components e2-NW, e2-E, 
and e2- W are marked with "+". The A3. 6 cm position of the UC HII region in 
W5 le2 Gau me et all 119931) is marked with " X ". 

Figure QJ the symbol "+" marks the positions of the compo- 
nents). The brightest source, W51e2-E, is located — 0.9" east 
of the UC HII regio n which is marked w ith " x " for the peak po- 
sition at 3 . 6 cm from lGaume et al.l (1 19931) . The secondary bright 
source, W51e2-W, coincides with the 3.6 cm peak of the UC 
HII region. The third component, W51e2-NW, is located about 
2" northwest of W51e2-E. In addition, these compact emission 
components appear to be surrounded by an amorphous halo in 
the continuum emission at 0.85 mm. The three bright contin- 
uum sources were also detected in the low-resolution contin- 
uum image at 1.3 mm (Figure |2j a)). Moreover, a weak contin- 
uum component, W51e2-N, is present — 2" north of W51e2-E. 

2.2. The VLA continuum images at 7 and 13 mm 

The observations of W51e2 at 7 and 13 mm were carried 
out with the VLA (see the summarized observing parameters 
in TableQJi. We made the calibrations by following the standard 
data-reduction procedure for VLA data with AIPSQ- Then, the 
calibrated visibilities were loaded into the Miriad environment. 



3 http.V/www. cfa. harvard, edii/sma/miriad 



4 http://www.aips.nrao.edu 
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FIG. 2. — (a): 1.3 mm continuum image of W51e2 observed with the SMA. Contours are ±5(7 X 2" (n = 0, 1, 2, 3, ... and a = 14.0 mJy beam -1 ). The FWHM 
beam (1.35" X 0.69", P.A.=— 87.4°) is shown at bottom right. The three components detected at 0.85 mm together with the additional component e2-N detected at 
1.3 mm are marked with "+". The position "+" of the submillimeter components and the peaks of the 3.6 cm continuum ("x") and H26a ("if") are also marked 
in the rest of the images, (b): 7 mm continuum image of W51e2 observed with the VLA. Contours are ±5ct X 2"/ 2 (n = 0, 1, 2, 3, ... and a = 1.12 mJy beam" 1 ). 
The FWHM beam (0.40" X 0.15", P.A.=80.3°) is shown at bottom-right, (c): 13mm continuum image of W51e2 observed with VLA. Contours are ±5cr X 2" 
(n = 0. 1, 2, 3, and a = 0.50 mJy beam" 1 ). The FWHM beam (0.23" X 0.10", P.A.=82.6°) is shown at bottom right, (d): Image of the H26a line emission 
integrated from 25 km s" 1 to 84km s~ ] from W51e2. Contours are ±5ct x 2" (n = 0, 1, 2, 3, ... and a = 0.146 Jy beam" 1 km s" 1 ). The FWHM beam (0.25") is 
shown at bottom right. All the RRL images are convolved to a common circular beam (0.25"). The star denotes the peak position of the H26a line and is marked 
in the rest of the RRL images, (e): Image of the H53c< line emission integrated from 22km s _1 to 101 km s~'. Contours are ±5cr X 2" (n = 0, 1, 2, 3, ... and 
a — 9.68 X 10~ 3 Jy beam -1 km s _1 ). (f): Image of the H66o line emission integrated from 24 km s~' to 120 km s~' . Contours are ±5cr x2"/ 2 (n = 0, 1,2, 3, ... and 
(j = 6.14 X 10~ 3 JybeairT 1 kms 4 ). 



The imaging and further analysis were carried out with Miriad. 

The dirty images were made using INVERT with robust weight- 
ing (robust=0) and were cleaned with the hybrid deconvolution 
algorithm. Figures|2jb) and|2jc) show the continuum images at 
7 and 13 mm, respectively. Strong continuum emission is de- 
tected from W51e2-W at both 7 and 1 3 mm in good agreemen t 
with previous observations at centimeters (Gaume et al. I H993h . 
In addition, a weak but significant (>10cr) emission feature at 
W51e2-E has been detected at only 7 mm. As shown in Figure 
|2jb), the weak 7 mm continuum source at W51e2-E has been 
clearly separated from the UC HII region, W51e2-W. 

2.3. Hydrogen recombination lines 

For the line cubes, the continuum emission was subtracted 
using the linear interpolation from the line-free channels with 
UVLIN. The H26a line (353.623 GHz) was included in the line 
cubes made from the SMA data at 0.85 mm. We also made a 
line image cube to cover the H30a line (23 1 .901 GHz) from the 
SMA data at 1.3 mm. Both the H26a and H30a line data were 
resampled to 1 km s -1 . The rms noises of 49mJy beam -1 and 
8 1 mJy beam" 1 in each of the channel images are inferred for 
the H26a and H30a line image cubes, respectively. Figure|2jd) 
shows the integrated line intensity image of the H26a line made 



with a 3(7 cutoff in the velocity range from 25 to 84km s . 
The H26a line emission from the UC HII region (W51e2-W) 
appears to be very compact and has not been resolved with the 
beam of 0.25", suggesting that the intrinsic source size of the 
hyper-compact HII core is < 0.06". The peak position of the 
integrated H26a line emission (stars in Figure |2]i appears to 
have a significant offset of ~0. 15" from the continuum peak at 
13mm (also see the Gaume's position in Figure|2td)). Unlike 
the other RRLs and continuum emission at centimeter wave- 
lengths, the line emission at the frequency of the H30a line 
(the image is not shown in this paper) shows an extended distri- 
bution covering the entire 3" region, suggesting that the H30a 
line might have been severely contaminated by one or more 
molecular lines in W51e2 (e.g., C3H7CN with a rest frequency 
of 23 1.9009 GHz). 

The line of H66a is broad and weak, which is barely covered 
by the VLA correlator band with which we had difficulty de- 
termining the line-free channels. Instead, the continuum levels 
were determined by averaging all channels excluding a few rel- 
atively strong line emission channels and were subtracted from 
the visibility data. The uncertainty is < 5% in the continuum 
levels due to the line contamination. The image cubes of the 
H53a and H66a lines were cleaned and convolved to a circular 
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Solid: H26o (hann=3) 
Dash-dots: H53a (x8.2, hann = 3) 
Dots: H66a (x15.B, hann = 3) 




15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 
Velocity (km/s) 

FIG. 3. — Profiles of hydrogen recombination lines, H26o, H53o, and 
H66«, integrated from W51e2-W in an area of 0.5" xO. 5" centered at the 
H26« line peak. All the line images have been convolved with the same cir- 
cular beam of 0.25". The spectra are multiplied with the ratios of the rest 
frequency of the H26a line to the RRL rest frequencies (vH26a/ v RRt)- The 
vertical dashed line denotes the systematic velocity of 53. 9± 1.1 km s which 
we determined from the previous measurements of seven different hot molec- 
ular lines (see Section 2.4). 



beam of 0.25", giving the rms noises of 2.4 and 1 .5 mJy beam" 1 
per channel, respectively. The integrated intensities of H53a 
(with a 3cr cutoff in each channel) and H66a (with a 2er cutoff in 
each channel) are shown in Figures [2je) and|2jf), respectively. 
Both H53a, and H66a lines are detected from only W51e2-W, 
the UC HII region. For the H26a, H53a, and H66a lines from 
W51e2-W, we integrated the line emission region with a size of 
0.5" around the H26a peak position. The profiles of the H26a, 
H53a and H66a lines are multiplied by the frequency ratios 
of v mba jv mba , v mba jv H ^ a and Vh26cl/vh(&c>, respectively, 
as shown in Figure [3] The measurements of the peak intensity 
Q peak), radial velocity (V LSR ), and FWHM line width (AV) for 
the H26a, H53a, and H66a lines are given in Table[3] 

2.4. Absorption ofHCN(4-3) and CO(2-l ) 

The systematic velocity of W51e2 vari es when it is mea- 
sured using differen t hot molecular lin es (Sollins et al. 2004; 
iRemiianet all 12004 IZhang etalj[l998t iRudolph et aiT Tl990). 
In order to minimize the possible effects of the different distri- 
butions with different molecular lines, we averaged these mea- 
surements of molecular lines H 13 CO + , S0 2 , SiO, CH 3 CN, and 
HCO + and obtained a mean value of 53.9±l.lkm s . We 
adopt this value as the systematic velocity for W51e2-E here- 
after. 

The lines of HCN(4-3) at v> = 354.505 GHz and CO(2-l) 
at vq = 230.538 GHz were covered in the SMA observations 
at 0.85 mm and 1.3 mm, respectively. The dirty images were 
made with natural weighting and deconvolved with the Hog- 
bom Clark Steer hybrid algorithm. The final rms noises of a = 
65 and 94mJy beam" 1 per channel (1km s" 1 in the channel 
width) were achieved for HCN(4-3) and CO(2-l), respectively. 
Figure 2] shows the spectra of HCN(4-3) toward W51e2-E and 
W51e2-W and the spectrum of CO(2-l) near the peak position 
of W51e2-E. A significant absorption of the HCN(4-3) line is 
detected against the submillimeter core, W51e2-E, and most of 
which is redshifted with respect to the systematic velocity of 
53.9±1.1 km s" 1 as indicated by the vertical dashed line in Fig- 
ure Uc). Note that the redshifted absorption of HCN observed 
with the high angular resolution against the compact continuum 



core corresponds to the cold gas in front of the dust core mov- 
ing toward it, i.e., the infall. The CO(2-l) spectrum shows a 
broader redshifted feature in absorption than that of HCN. Due 
to the poor angular resolution in the CO(2-l) observations, the 
observed broad profile in absorption is caused not only by the 
infall but also by the outflow. In particular, the high redshifted 
absorption arises very likely from the outflow gas. In compari- 
son with that of W5 Ie2-E, the spectrum toward W5 Ie2-W only 
shows a weak spectral feature (Figure|4|d)) which can be char- 
acterized as an inverse P Cygni profile, suggesting that only 
little molecular gas (if present) falls onto the UC HII region. 

To show the spatial distribution of the absorption gas, we car- 
ried out a moment analysis. The calculation of the zeroth mo- 
ment corresponds to the integrated intensity over the velocity. 
In order to avoid cancellation between emission and absorption 
spectral features in the same beam area, we separate the absorp- 
tion and emission in the moment calculation. We found that 
both the blueshifted and redshifted emission of the HCN(4-3) 
line (the images are not shown in this paper) is extended from 
northwest to southeast across W51e2-E, which is, in general, 
consistent with the molecular outflo w direction as suggeste d 
based on the CO(2-l) line observation ( Ket o & Klaassenl20 08). 
The HCN(4-3) absorption, integrated from 44 to 62 km s with 
a 5a cutoff in each channel (FigureUa)), is found to be mainly 
concentrated on W51e2-E. From the lower-resolution observa- 
tion, the distribution of the CO(2-l) absorption line integrated 
from 43 to 71 km s" 1 (a 5a cutoff in each channel, see Figure 
|4|b)), also shows the absorption peaks at W51e2-E instead of 
W51e2-W. We also note that the CO(2-l) absorption shows a 
relatively extended feature north of W51e2-E, further indicat- 
ing that the absorption of CO(2- 1 ) might indeed be significantly 
contaminated by the outflow gas. 

Along with the absorption spectra, the distribution of the ab- 
sorbing gas in W51e2 evidently demonstrates that the submil- 
limeter (dust) core W51e2-E, instead of W51e2-W (the UC HII 
region), is the center of accretion for the majority of the high- 
density gas from this molecular core. 

2.5. Spectral Energy Distributions 

In lack of high angular resolution observations at millime- 
ter/submillimeter wavelengths, the flux densities from the en- 
tire region of W51e2 were considered in the previous analyses 
of the SED (IRudolph et al.ll 1 990t [Sollins et alj2004tlKeto et al.1 
2008). With the high-resolution observations at the wavelengths 
from radio to submillimeter, we are now able to spatially sepa- 
rate individual emission com ponents in our SED analysis. 

Addi ng the data at 3.6 c m dGaume et a n ri993h and at 2 and 
6 cm ( Rudolp het al.l Il990h along with the continuum data at 
the four wavelengths (13, 7, 1.3, and 0.8 mm) discussed in this 
paper, we have seven measurements in flux densities covering 
nearly 2 order of magnitude in frequency. For 13, 7, and 0.8 
mm data, the flux density measurements were made by fitting 
a Gaussian to the individual emission components in the im- 
ages with the same size of the convolved FWHM beam (0.4"). 
For the 1.3 mm data, we determined the flux density by fit- 
ting multiple Gaussian components to the emission features in 
a relatively lower resolution (1.1"), which results in relatively 
large uncertainties. For the emission cores W51e2-W, W51e2- 
E, W5 le2-NW, and W5 Ie2-N, the flux densities determined are 
given in Table [2] The relevant continuum spectra for the four 
emission components are shown in Figure[5] 
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FIG. 4. — (a): Image of the HCN(4-3) absorption line integrated from 44 to 62 km s . Dashed contours are ±5<r X 2" (n = 0, 1, 2, 3, ... and 
a = 0.113 Jy beam -1 km s~'), overlaid on the gray-scaled continuum image at 0.85mm. The FWHM beam 0.33" X 0.24" (sup=0) is shown at bottom right. 
(b): Image of the CO(2-l) absorption line integrated from 43 to 71 km s -1 . White dashed contours are ±5o- X 2" (n = 0, 1, 2, 3, ... and cr = 0.771 Jy beanr' km s" 1 ), 
overlaid on the gray-scaled continuum image at 0.85 mm. The FWHM beam 1.35" X 0.73" (sup=0) is shown at bottom right, fc and d): The spectral profiles of 
the HCN(4-3) absorption line toward W51e2-E and W51e2-W, respectively. The vertical dashed line denotes the systematic velocity of 53.9±1.1 km s . (e): The 
spectral profile of the CO(2-l) absorption line toward W51e2-E, with a vertical dashed line for the systematic velocity. 



3. NATURE OF THE W51E2 COMPLEX 

Four distinct components have been identified in the 3" re- 
gion of W5 le2 from the high-resolution images at wavelengths 
from centimeter to sub millimeter (Figure Q] and [2]). The nature 
of these components is discussed as follows. 

3.1. W51e2-W: UCHII region 

W5 Ie2-W is the only core detected at centimeter wavelengths, 
and its continuum emission is dominated by the free-free emis- 
sion from the ionized core and traces the thermal HII region 
around the central star. Figure |2jb)-(f) show the UC HII re- 
gion, W51e2-W, including an unresolved core (a possible ion- 
ized disk) and a northeast and southwest (NE-SW) extension 
as a possible outflow. The NE and SW extensions are best to 
be seen at 7 and 13 mm in Figure|2jb) and (c), which is in good 
agreement w i th the previous observation s at 3.6 and 1.3 cm by 
iGaume et ail dl993l) . iGaume et all d 19931) found the spectral in- 
dices a of about 2 (optically thick) and 0.4 for the NE and SW 
extensions of W51e2-W, respectively, and suggested that there 
is a one-side (SW) collimated ionized outflow from the core. 

The H26ct line is an excellent tracer for hyper-compact ion- 
ized cores (star in Figure |2fd)), providing a good diagnosis for 
the presence of an ionized disk. On the basis of the SMA 
observations with an angular resolution of 0.25", the hyper- 
compact ionized core which peaked at 59.1 km s -1 in W51e2- 
W has not been resolved, giving a limit on the intrinsic size of 
9 S < 0.06" and linear size < 310AU from a Gaussian fitting. 
The elongation shown in the H53a image (Figure|2je)) agrees 



with the radio continuum images observe d at 1.3 cm (F igure 
12c) and IGaume et"aT]|1993l) . 3.6 cm dGaume et al.lll993h . and 
7 mm (Figure |2b)), suggesting that the NE-SW extension ob- 
served in the H53a line corresponds to the expansion of ion- 
ized gas (or outflow). The H66a line emission also has a sim- 
ilar elongation (Figure [2f)), with peak intensity near Gaume's 
position in the NE and an extension structure in the SW. We 
also found a significant velocity gradient of the H53ct line in 
the W5 Ie2-W region, in agreement with what was observed by 
iKeto & Klaasserj d2008). redshifted in the SW and blueshifted 
in the NE. The broad wings of the recombination lines shown 
in Figure [3] (also in Table [3) corresponding to extended, opti- 
cally thin emission indicate that it is hard to explain them as an 
ionized disk. The velocity gradient observed in the H53a line 
is more likely produced by an ionized outflow rather than an 
ionized accretion disk. 

We fitted the SED of W51e2-W with a free-free emission 
model 

S„ = n s B„(T e )(l-e- T <) (Jy), (1) 
where Si, is the flux density at frequency v, tt s is the solid angle 
of the source, Bv(T e ) is the Planck function, and T e is the elec- 
tron temperature. The contin uum optical depth is expressed by 
dMezger & Henderson! 19671) as 

fT e \~ 135 f v \-2.i / EM \ 

The emission measure (EM) can be written as 

EM = A£L/v-, (3) 
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Table 2 

Core fluxes of W5 1e2 at Centimeter to Submillimeter bands 



W51e2-W W51e2-E W51e2-NW W51e2-N References and Beamsize 6 



parameters 



R.A.(J2000) 
Decl.(J2000) 
Deconvolved size (" 

So.85mra(Jy) 

Sl.3mm(Jy) 

S7mm(Jy) 

Sl.3cm(Jy) 

S 2 cm(Jy) 

S3.6cm(Jy) 

SficmCJy) 



19:23:43.90 19:23:43.96 19:23:43.87 19:23:43.97 This paper 
+ 14:30:34.62 +14:30:34.56 +14:30:35.97 +14:30:36.57 This paper 



0.45 
0.35±0.07 

0.62±0.09 
0.43±0.03 
0.21±0.05 
0.067±0.004 
0.025±0.003 



0.71 
3.30±0.20 
2.15±0.12 
0.04±0.01 

< 0.004 

< 0.001 



0.38 
0.81±0.34 
0.62±0.12 

< 0.0052 

< 0.004 

< 0.001 



1.24 

0.73±0.08 

< 0.0052 

< 0.004 

< 0.001 



This paper 
This paper: 
This paper: 
This paper: 
This paper: 



9 = 0.4" 
9 = 1.1" 
6 = 0.4" 

6 = 0.4": Gaumeetal. (1993) : upper limit 
Rudolph et al. (1990) : 6 = 0.4" 

Gaume et al. (1993) : 6 = 0.21": Gaume et al. (1993) : upper limit 
Rudolph et al. (1990): 6 = 0.6" 



*The sizes of e2-W, e2-E, and e2-NW are derived from the 0.85 mm image, and e2-N is from the 1.3 mm image. 



where N e is the mean electron density, L is the path length, 
and fv is the volume filling factor. The dimensionless f actor 
a(y, T e ) is the order of unity (Mezger & Henderson 1967). For 
a homogeneous ionized core, the model can be expressed with 
two free parameters, namely, T e and EM, for the first approxi- 
mation described above. 

The solid curve in Figure [5j a) shows the best fit to the data, 
suggesting a turnover frequency of v$ ~ 27 GHz. The derived 
mean values for the physical parameters for the overall HII 
region, 7>4900±320K and E M=(1.3±0.2)x 10 9 p c cm" 6 , are 
similar to the previous result in Zhan g"etai] d!998l) . Thus, the 
mean electron density of N e ^3xl0 5 cm" 3 is inferred on the 
assumption of Lf v ^0.5"x5.1 kpc. 

Both the Lyman continu um flux (Nl) and the excitation pa- 
ramet er (£/) can be derived dRubinll96 8: Panagi all973tlMatsakis 
119761) as follows: 



Assuming N(He) /N(H) = 0.096 dMehringerl 19941) and a(v, T e ) = 
1 , we have T* = 12, 000 ± 2 , 000 K which appears to be consid- 
erably higher than the mean electron temperature of the overall 
HII region. The EM of the hyper-compact HII core can be as- 
sessed as 

3/2 / n N -2 



EM = 7.1pc cm -0 [ ^ 



Ay 



A 



km s 



(7) 



where T e is the electron temperature, S is the intrinsic size of 
the HII region, Sl is the peak line intensity in Jy, AV is the 
FWHM line width in km s -1 , and A is the observing wave- 



^— jlength in millimeters. Based on the measurements of the H26a 
line from the hyper-compact ionized core together 



3gether with the as- 



jVz.>7.5xl0 46 s-' \^-\ ( — V 



Jy J \kpc J 



(gHz) 



o.i 



\ 10 4 K J 



(4) 
(5) 



where D is the distance to the source. We assumed a(v, T e ) ~ 1. 
Based on our model fitting, the continuum emission from e2-W 
becomes optically thin at a frequency >80GHz (r <0.1), and 
we obtained A^=3.0xl0 48 s" 1 and t/=37.6pc cm -2 . The in- 
ferred Lyman continuum flux requires a massive star equivalent 
to a zero-age mean-sequen ce star of type Q8 located inside the 
UC HII region W51e2-W (Panagia fl 19731). wh ich is consistent 
with the result of 07.5 in lRudolph et al.ldl990h . Alternatively, a 
cluster of B-type stars can also be responsible for the ionization. 
The total ionized mass within 0.5" of W51e2-W is ~0.02M W , 
which agrees with the value derived by Ga ume et alj (Il993l) . 

For the hyper-compact core (9 S < 0.06") as observed with the 
H26a line, the electron temperature can be estimated from the 
H26a line and the continuum flux density at 354 GHz (0.85 mm) 
on the assumption of optically thin, LTE condition: 



y * 




So 



85 r 



Sh26o 



i , N(He) 
. N(H) 



(6) 



sumption of T e w T*, EM > 7 x 10 10 pc cm -6 is found. The 
corresponding lower limit of the volume electron density is 
N e > 7 x 10 6 cirT 3 assuming Lf v < 0.06" x 5.1 kpc. The in- 
ferred high electron temperature and high electron density sug- 
gest that the H26a line arises from a hot, very compact re- 
gion (< 0.06") which is probably close to the central ionizing 
star(s). With the nature of optically thin and sensitive to the 
high-density ionized gas, the H26a line is an excellent tracer 
of the hot, high-density ionized region surrounding the ioniz- 
ing source. In comparison with the mean electron tempera- 
ture (7^,=4900 K) inferred from the free-free emission in a large 
area (0.5"), the high temperature of 7^=12,000 K derived from 
the H26a in a compact area (< 0.06") suggests that a temper- 
ature gradient is present along the radius of the UC HII region 
W51e2-W. 

Observations of RRLs in a wide range at wavelengths from 
centimeters to submillimeters also offer a means to further ex- 
plore physical conditions of W51e2-W. Table[3]summarizes the 
reliable measurements of three RRLs from W51e2, namely, 
H26a, H53a, and H66a. We plotted the spectral profiles of 
the three RRLs together, and each of the profiles has been mul- 
tiplied by a factor of V2ba/v„ a , the ratio of the H26a line fre- 
quency to that of RRLs (Figure |3)- In the case of optically 
thin and no-pressure broadening, the three line profiles should 
match each other. Comparing H53a with H26a, we find that 
the peak intensities of the two lines agree with each other, which 
suggests that, in general, both the H53a and H26a lines are un- 
der an optically thin, LTE condition. Several velocity peaks in 
the H26a line indicate multiple kinematic components in the 
hyper-compact HII core. The H53a line is characterized by a 
relatively smooth profile with large velocity wings. The large 
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FIG. 5. — The SEP for W51e2 cores, (a ): SE P fitting for W51e2- W with a free— free emission model. Black points denote the measurements from this work, and 
the asterisks are from lRudolph et alj fT990[) and lGaume eTal] ft993l). All flux values are lis ted in Table[2] (b): The SEP fitting for the W51e2-E core with a thermal 
dust emission model. The arrows denote the upper limits from iRudolph et all ft990l) and Gaume et al. 1 1993). (c): The SEP fitting for the W51e2-NW core. The 
upper limit of the flux density at 43 GHz (7 mm) is from our data (3<r of the 7 mm image), (d): The SEP fitting for the W51e2-N core. 



Table 3 

Parameters of RRLs derived from our data* 



velocity wings observed in the H53a line can be explained by 
the large velocity gra dient in the lower-d ensity electron gas of 
the ionized outflow (Gaume et al. 1 119931) . However, the line 
intensity of H66a appears to be significantly weaker than the 
other two high-frequency RRLs. The peak of the modified line 
profile of H66a is a factor of ~ 5 less than that of the other 
two. Considering the fact that the frequency of H66a is be- 
low the turnover frequency of i/q ks 27 GHz, we found that the 
H66a line from the ionized gas is obscured severely due to the 
self-absorption process in the hyper-compact HII core. From 
Equation (2), a mean optical depth of t c = 1 .6 at VH66a = 22.36 
GHz is found. The exponential attenuation of exp(-1.6) « 0.2 
suggests that the H66a line is attenuated mainly due to the self- 
absorption in the ionized core. In addition, the pressure broad- 
ening effect which weakens the lower frequency lines needs to 
be assessed. 

The total line width of RRLs (Ay) can be expressed by the 
Doppler b roadening (AVp) and pres sure broad ening (AV P ). Ac- 
cordin g to lBrocklehurst & Leemanld 19711) and lGordon & Sorochenkot vhere AV, is the turbulent broadening of the gas including broad- 





H26a 


H53q 


H66a 


Rest Frequency (GHz) 


353.623 


42.952 


22.364 


Measurements 








Ipeak (Jy) 


1.29±0.02 


0.137±0.001 


0.015±0.003 


Vlsr (km s- [ ) 


59.1±0.2 


60.4±0.2 


62.1±1.3 


AV (kms -1 ) 


23.1±0.5 


32.0±0.3 


35.1± 8.3 


Derived line broadenings 








AVp (kms" 1 ) 


0.01 


1.3 


6.5 


AVr (km s" 1 ) 


15.0 


15.0 


15.0 


AV, (kms -1 ) 


17.5 


28.3 


31.1 


AV D (km s" 1 ) 


23.1 


32.0 


34.5 



*A11 line profiles are integrated from the 0.5" xO. 5" region around 
the peak position of H26a from maps of the same beamsize of 0.25" . 



(2002). we have 



AV D = km s" 1 4 0.0458 x 



AV P = 3.14 x 10" 14 kms 



7; 

K 

-1 -4.4 



/ AV, 



\km s 



mm 



N e 



cm 



-o.i 



Ay 



Ay^+Ay2, 



ening due to outflows and disk rotations and n is the princi- 
pal quantum number of the transition. We calculated broaden- 

(g) ings on the mean electron temperature (T e «4900K) and the 
mean electron density (N e «3 x 10 5 cm" 3 ) of the UC HII region. 
The thermal broadening (AV r ) is 15 km s" 1 for 7>4900K. The 
pressure broadening is a strong function of the principal quan- 
tum number, i.e., AVp increases drastically toward the low- 
frequency lines. For the high-frequency line (H26a), the pres- 

* ' sure broadening (AVp = 0.01 km s" 1 ) can be ignored, while for 
the H66a line, the pressure broadening of AVp=l km s be- 
(10) comes significant. The pressure broadening is Ayp=76km s" 1 
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for the H92a line, which can actually wash out the profile of the 
line emission. Thus, in addition to the opacity effect, the pres- 
sure broadening may also make a considerable contribution to 
diminis h the H92a fine w hich was not detected in the observa- 
tions of Mehr ingeri d 1994b . The turbulent and/or the dynamical 
motions of the ionized gas contribute a total of AV, «18 km s -1 
in the line broadening for the H26a line and AV, «28 and 
31km s" 1 for the H53a and H66a lines, respectively. The 
H26a line traces the high-density ionized gas in the hyper- 
compact core or the ionized disk (if it exists) and is less sen- 
sitive to the lower-density outflow which produces the broad 
wings in the line profiles. The turbulent and/or dynamical broad- 
ening becomes dominant in the H53a and H66a lines, and 
these lower frequency RRLs (H53a and H66a) better trace the 
lower-density electron gas in the ionized outflow. 

In short, W51e2-W consists of a hyper-compact HII core (< 
310 AU) with an emission measure of EM > 7 x 10 10 pc cm" 6 
and an ionized outflow. An early-type star equivalent to an 08 
star (or a cluster of B-type stars) is postulated to have formed 
within the hyper-compact HII core. The broadening of the line 
profiles is dominated by the Doppler broadening including both 
thermal and turbulence/dynamical motions. Pressure broaden- 
ing becomes significant only for the H66a line and the lines 
with a larger principal quantum number. 

3.2. W51e2-E: A massive proto-stellar core 

W5 Ie2-E is located - 0.9" east of the UC HII region W5 le2- 
W. It is the brightest source in the 0.85 mm continuum image 
(Figure [T]). The flux density drops drastically at 7 mm (Figure 
|2jb)). It cannot be detected at longer wavelengths. The SED of 
the continuum emission from W51e2-E (Figure |2b)) appears 
to arise from a dust core associated with proto-stars. 

With the Rayleigh-Jeans approximation, the continuum flux 
density fro m a homogeneous, isotherm al dust core can be de- 
scribed as dLaunhardt & Henn ing 1997) 




where kq is the dust opacity at frequency i/q, f3 is the power- 
law index of the dust emissivity, is the mass of dust, 7^ is 
the dust temperature in K, and D is the distance to the source. 
We adopt the dust opacity value of k(X = 1 3mm) = 0-8 cm 2 g -1 
dOssenkopf & Henning|l994l:lLaunhardt & Henningll 19971) and 
have three free parameters for the dust core (Mj, Tj, and j3). 

The solid line in Figure |5fb) shows the best least-squares fit- 
ting to the observed data. The index of /3 = 0.26 ± 0.08 de- 
rived from our fitting indicates that emission from the e2-E core 
migh t be dominated by large grains of du st ( Miyake & Nakagawa 
I1993L ICheneTaT]|1995l:lAndrews& Williamsll2007h . In addi- 
tion, because of a lack of flux density measurements at higher 
frequencies or shorter wavelengths, the dust mass and tempera- 
ture are degenerate i n our best fitting, i.e., M^Td = 142±9M K. 
Rem ijan et alJ d2004l) derived the kinetic temperature in W5 le2 
to be 153 ±21 Kfromthe high-resolution observations of CH3CN, 
which is close to the dust temperature 7^=100 K used in Zha ng et ail 
d!998l) . Adopting that T d =100K and assuming that the ratio of 
H2 gas to dust is 100, we found that the total mass of ^14OM 
is in the W51e2-E core. Emission from the W51e2-E region 
has been slightly resolved, so that the dust core of W51e2-E 
might host a number of proto-stellar cores which accrete the 



surrounding molecular gas as indicated by both the HCN and 
CO absorption lines (see Section 2.4). This proposed scenario 
for W5 Ie2-E is also consistent with the hourglass-lik e magnetic 
fields inferred from polarization measurements of Tang et alJ 
(2009) for the W51e2 region. The configuration center of B 
vectors coincides with the peak position of the W51e2-E dust 
core. 

From a Gaussian fitting to the absorption spectrum of the 
HCN(4-3) toward the center of W51e2-E (FigureUc)), we de- 
termined the peak line intensity All = —1.3 ±0.1 Jy beam" 1 , the 
line center velocity Vhcn =56.4 ± 0.2 km s" 1 and the line width 
of AVhcn=9-5 ± 0.4km s" 1 . The infall velocity of the gas is 
the offset between the center velocity and the systematic veloc- 
ity, i.e., Vi„ = VHCN-V S y S « 2.5 km s -1 . The peak continuum 
intensity of W51e2-E at 0.85 mm is I c = 1.2 ± 0.01 Jy beam -1 
(corresponding to a brightness temperature of 144 ± 1 K). The 
ratio of -All/ Ic ~ 1 suggests that the absorption li ne is satu- 
rated. Thus, a lower limit on the optical depth (see lOin et al] 
2008) is thcn > 3. Assuming that the excitation temperature 
of HCN(4-3) equals the dust temperature, 100K, we obtained a 
lower limit of HCN colu mn density of 7.9 x 10 15 cm" 2 . Taking 
[H 2 ]/[HCN] ~ 0.5 x 10 8 dlrvine etal.fl 987) and the size of the 
infall region as 0.75" (~ 4000 AU, see Figure Ufa)), we found 
that the hydrogen volume density Nh 2 ~ 6.9 x 10 6 cm" 3 . The 
infall rate of the gas can be estimated by 

dM/d? = 2.1 x lO^Moyr" 1 ( J*-) ( JL) ( -^ T ) x 

\arcsecy \kpc/ ykms -1 / 

(^H) - (12) 
V 10 6 cm" 3 J ' 

where 20 m = 0.75" is the diameter of the infalling region, D 
is the distance to the source, V,„ is the infall velocity, and the 
molecular mass ratio mjmu^ = 1.36 is assumed. With the de- 
rived parameters a lower limit of infall rate 1 .3 x 10~ 3 M© yr" 1 
is inferred, suggesting that W51e2-E is the accretion center of 
the W51e2 complex. 

3.3. W51e2-NWandW51e2-N 

W5 le2-NW was clearly detected bv lTang et all d2009l) in the 
continuum emission at 0.87 mm. They also found a signifi cant 
concentration of B fields at the e2-NW. lGenzel et al.1 dl98ll) and 
Imai et all d2002l) detected bright and compact H2O masers near 
this source. We also detected local enhanced continuum emis- 
sion at 0.85 and 1.3 mm. No significant continuum emission 
has been detected at 7 and 13 mm (Figur^U and |2), and there 
are no significant absorption lines of HCN(4-3) and CO(2-l) in 
this region. We fitted the SED of continuum emission with a 
thermal dust model as shown in Figure 0c). A lower limit of 
f3 > 0.35 is derived from the fitting, which leads to an upper 
limit of MjTj < 39 M Q K. Assuming the dust temperature in 
this region to be ~100K, we inferred an upper limit of 40 M© 
in the total gas mass if the ratio of H2 gas to dust is equal to 
100. W51e2-NW appears to be a massive core at a very early 
phase of star formation. 

Located ~ 2" north from W51e2-E, W51e2-N is detected 
in continuum emission from our lower-resolution (1.4" x 0.7") 
image at 1.3 mm (Figure |2j a)). This source appears to be re- 
solved out with a high angular resolution (0.3" x 0.2") obser- 
vation at 0.85 mm. No significant continuum emission was de- 
tected at the longer wavelengths. The SED consisting of the 
flux density at 1.3 mm and the upper limits at 7, 13, and 36 mm 
is also fitted with a thermal dust model. From the fitted val- 
ues of the parameters (3 > 0.75 and M^Td < 68 M Q K, we found 
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an upper limit of M d < 70 M if T d ~100K and the H 2 gas to 
dust ratio is 100. No maser activities have been detected from 
W51e2-N, suggesting that W51e2-N is probably a primordial 
molecular clump in W51e2. 

3.4. A propagating scenario of star formation in W51e2 

Apparently, in the massive molecular core W51e2, the gravi- 
tational collapse occurred first at W5 Ie2-W where an 08 star or 
a cluster of B-type stars were formed. Based on the absorption 
spectrum of the HCN line, there appears to be little molecular 
gas present in the immediate environs. The accretion appears 
to be paused by the intensive radiation pressure from the cen- 
tral ionized star(s). Most of the surrounding gas (^0.02 M ) 
has been ionized. The thermal pressure in the UC HII region 
drives an ionized outflow from the ionized core. The submil- 
limeter observations of dust emission and molecular line ab- 
sorption show that the major accretion now has been re-directed 
to W51e2-E, 0.9" (~4600AU) east from W51e2-W. W51e2-E 
becomes the new dominant gravitational center accreting mass 
from the surroundings with a rate of 10" 3 M Q yr" 1 . Star for- 
mation activities take place in this massive core, as shown by 
maser activities, outflow, and organized B-field structure. Since 
no free-free emission and no RRLs have been detected, the 
massive core W51e2-E likely hosts one or more massive proto- 
stars. The offset of the radial velocities between the two cores, 
W51e2-E and W51e2-W, is SV = 59.1 -53.9 = 5.2km s" 1 . If 
W5 Ie2-W and W5 Ie2-E are gravitationally bound and W5 le2- 
W is circularly orbiting around W51e2-E, W51e2-E has a dy- 
namical mass of > 140 M , in good agreement with the mass 
derived from the SED analysis. 

A bipolar outflow from W51e2-E has been suggested from 
the observations of molecular lines. The impact of the outflow 
on the medium in its path may trigger the collapse of the sub- 
core W51e2-NW and induce further star formation activities 
there. W51e2-N represents a gas clump with a considerable 
amount of mass for the star formation. This speculative sce- 
nario appears to reasonably explain what we have observed in 
the W51e2 complex. 

4. SUMMARY AND CONCLUSIONS 

We have presented high-resolution images of the W51e2 re- 
gion at 0.85, 1.3,7, and 1 3 mm for continuum and hydrogen re- 
combination lines (H26a, H53a, and H66a) and the molecular 
lines (HCN(4-3) and CO(2-l)). The W51e2 complex has been 
resolved into four distinct components, W51e2-W, W51e2-E, 
W51e2-NW, and W51e2-N. We have carried out a comprehen- 
sive analysis of the continuum SED, the RRLs, and the molecu- 
lar absorption lines and found that the four cores are at different 
phases of massive star formation. 

1. W51e2-W, associated with the UC HII region, is the only 
source from which the H26a, H53a, and H66a lines have been 
detected. The Lyman continuum flux inferred from the SED 
analysis suggests that a massive 08 star or a cluster of B-type 
stars have been formed within the HII region. The unresolved 
H26a line emission region suggests the presence of a hot (T* = 
12, 000 ±2, 000 K), hyper-compact ionized core with a veloc- 
ity of 59.1 ±0.2 km s" 1 , a linear size of < 310 AU, and a large 
emission measure EM> 7 x 10 10 pc cm" 6 (or a high density 
N e > 7 x 10 6 cm -3 ). The H53a and H66a line images show 
an SW elongation from the hyper-compact ionized core, sug- 
gesting either an expansion or an outflow of the ionized gas 
with relatively lower density. No significant detection of the 



HCN molecular line in absorption against the compact HII re- 
gion indicates that the significant accretion onto this core has 
been stopped. The line ratios between the H26a, H53a, and 
H66a show that both the H26a and H53a lines from W51e2- 
W are in an optically thin, LTE condition. The emission of 
the H66a line from the hyper-compact core appears to be at- 
tenuated mainly due to the self-absorption in this HII region. 
The line profile of the H26a appears to be dominated by the 
Doppler broadening due to the thermal motions of the hot elec- 
trons while the line broadening in both the H53a and H66a is 
dominated by the Doppler effect due to the dynamical motions 
of the ionized outflow. The pressure broadening in both the 
H26a and H53a lines is negligible and becomes significant for 
the H66a line corresponding to AVp <~ 7 km s" 1 . 

2. W51e2-E is a massive (~14OM ) dust core as suggested 
based on the SED analysis. No hydrogen recombination lines 
and no radio continuum emission (A > 1 cm) have been detected 
from this region, which suggests that W51e2-E is the major 
core hosting massive proto-stellar objects. Both the absorptions 
of HCN(4-3) and CO(2-l) are redshifted with respect to the 
systematic velocity, 53.9 ± 1.1km s" 1 , indicating that a large 
amount of molecular gas moves toward the massive core. From 
the ratios of the absorption line to continuum, a large infall rate 
of > 1.3 x lO~ 3 M yr _1 is inferred from the HCN absorption 
line, corresponding to the accretion radius of 2000 AU. 

3. W51e2-NW appears to be a massive (< 40 M ) core at 
an earlier phase of star formation. Given its location in the 
projected path of the molecular outflow from W5 Ie2-E, the star 
formation activities of W51e2-NW are likely to be triggered by 
the outflow impact on the medium in the W51e2 complex. 

4. W5 Ie2-N appears to be a relatively less dense but massive 
clump (< 70 M ) from which no obvious evidence has been 
found for star formation activities. 
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